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Abstract. The formation of globular clusters is still an unsolved prob- 
lem. Though most scenarios assume a massive molecular cloud as the 
progenitor, it is unclear, how the cloud is transformed into a star cluster. 
Here a scheme of supernova (SN) induced cluster formation is investi- 
gated. In this scenario the expanding SN shell accumulates the mass 
of the cloud. This is accompanied by fragmentation resulting in star 
formation in the shell. If this stellar shell expands sufficiently slow, its 
self-gravity leads to a recollapsing shell, by this forming one or several 
stellar clusters. 

In this paper N-body simulations of collapsing shells moving on cir- 
cular orbits in a galactic potential are presented. It is shown that typical 
shells (10 5 Mq, 30 pc) evolve to twin clusters in the galactocentric dis- 
tance range between 3 and 11 kpc. Their masses show a strong radial 
trend: on orbits inside 5 kpc both clusters have almost equal mass. Out- 
side 5 kpc the more massive twin cluster contains about 55% of the shell's 
mass, whereas the mass of the smaller decreases linearily to 15% at 11 
kpc. Outside 11 kpc the collapsing shells end up in a single cluster. Inside 
3 kpc the shells are tidally disrupted and only fragments substantially less 
massive than the initial shell survive. 

1. Introduction 

Several scenarios for the formation of globular cluster have been suggested, e.g. 
the collapse of giant molecular clouds (GMC) or the collision of molecular clouds 
(e.g. Fall & Rees 1985, Murray & Lin 1990, Fujimoto & Kumai 1997). A com- 
mon characteristic of all these scenarios is, that the clusters are formed from 
smooth gaseous distributions which are transformed into stars. However, this 
assumption requires short formation timescales and unusually high star forma- 
tion efficiencies in order to end up with a gravitationally bound cluster. An 
alternative model suggested by Brown et al. (1991) can overcome these difficuli- 
ties: their scenario starts with an OB-association exploding near the center of 
a molecular cloud. The expanding shell sweeps up the cloud material and in a 
later stage the expansion can be decelerated or stopped by both, the accumu- 
lated mass and the external pressure of the ambient hot gas. The shell itself 
undergoes fragmentation and, finally, forms stars. If these stars are gravitation- 
ally bound, they will recollapse, by this forming a globular cluster. 

Though a discrimination between the scenarios by means of hydro dynamical 
simulations (starting from first principles) is far out of reach at the moment, 
one can study different evolutionary stages. E.g. Theis (2000) compared in a 
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Figure 1. Temporal evolution of the masses of the twin clusters 
(signs) and the mass of unbound particles (solid line). Shown are re- 
sults for circular orbits at a galactocentric distance of 10 kpc (left) and 
5 kpc (right). The twins in the right panel collide at t ~ 28.6 and 
merge at t « 33.4. The total mass in each simulation is 10 5 M . The 
time unit is 7.7 Myr. 



series of N-body simulations the collapse of thin stellar shells and homogeneous 
spheres in a galactic tidal field. These calculations were performed for circular 
and eccentric orbits, but with a constant apogalacticon of 5 kpc. It was found 
that collapsing shells preferably end in multiple systems, mainly twins, whereas 
homogeneous spheres either form single clusters or become completely disrupted. 
The twins might survive for - at least - several galactic revolutions, and some 
of them show kinematical peculiarities like counter-rotating cores. 

In this paper I will address the question, how the formation of twins depends 
on the galactocentric distance, i.e. the strength of the galactic tidal field. 



2. Numerical Models 

The simulations start with a thin, spherical shell of a mass M to t = 10 5 M Q , an 
outer radius of 30 pc and a thickness of 3 pc. The shell is initially at rest in the 
sense that there is no overall expansion or contraction of the shell with respect 
to its center. The potential of the Galaxy is modelled by an isothermal halo with 
a circular speed of 220 kms -1 . In all simulations the initial velocity of the shell 
corresponds to a circular orbit. The calculations are performed with N = 10 4 
particles using a GRAPE3 board. 

Fig. 1 shows the temporal evolution of the masses of the two largest clusters 
for galacatocentric radii of 5 and 10 kpc. At 10 kpc two clusters with a mass 
ratio of 5:2 are formed after 20 Myr. Both clusters survive until the end of the 
simulation (~ 400 Myr). 20% of the stars initially residing in the shell became 
unbound. A slightly different behaviour is seen for the model starting at 5 kpc: 
Again twins are formed, but they have almost equal mass. After 210 Myr they 
collide and, finally, they merge. Both events are clearly reflected in step-like 
increases of the number of unbound particles which amounts finally to 40% of 
the mass of the shell (Fig. 1, right). 
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Figure 2. Cluster mass vs. galactocentric distance. Shown are the 
masses of the five most massive clusters of each simulation 18.6 Myr 
after the collapse. The distribution is unaffected by later merging of 
clusters. The clusters are sorted by mass (left: most massive clusters 
(box, plus), right: smaller clusters (triangle, filled circle, Y)). 



The masses of the formed clusters exhibit a clear radial trend (Fig. 2): The 
mass of the largest cluster increases almost linearily from 2 to 5.5 kpc reaching a 
plateau of about 55% of the total mass. Beyond 11 kpc only one massive cluster 
is formed. About 40% of the stars are not bound in clusters. Inside 3 kpc, a 
large set of small clusters is formed instead of a dominating pair of clusters. 
The mass ratio q between the two most massive clusters also shows two regimes: 
Below 5 kpc the masses are almost identical, whereas outside 5 kpc q increases 
up to ~ 4 : 1 at 11 kpc. In addition to the large clusters (> 10 M ), typically 
several small (< 10 3 M©), gravitational^/ bound clusters are formed. Contrary 
to the large clusters, their mass shows no trends with galactocentric distance. 

The simulations demonstrate that twin formation is expected over a large 
radial range. About 1/3 of these twins merge within 400 Myr after their forma- 
tion. The surviving twins are characterized by large spatial separations which 
makes them less likely to undergo a subsequent merger. Therefore, twin globu- 
lars might still exist in the Milky Way, but they can be unidentified as twins due 
to their large separation. Since they should share common orbits and metallici- 
ties, they could be found by proper determinations of globular cluster orbits. 
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